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We provide description of rapidity spectra of particles produced in pp collisions using anomalous diffusion ap- 
proach to account for their non-equilibrium character. In particular, we exhibit connection between multiproduc- 
tion processes and anomalous diffusion described through the nonlinear Focker-Planck equation with nonlinearity 
given by the nonextensivity parameter q describing the underlying Tsallis g-statistics and demonstrate how it 
leads to the Feynman scaling violation in these collisions. The q parameter obtained this way turns out to be 
closely connected to parameter l/k converting the original poissonian multiplicity distribution to its observed 
Negative Binomial form. The inelasticity of reaction has been also calculated and found to slightly decrease with 
the increasing energy of reaction. 



1. INTRODUCTION 

In description of multiparticle production pro- 
cesses one often uses statistical methods and con- 
cepts which follow the classical Boltzmann-Gibbs 
(BG) approach. However, it was demonstrated 
recently that to account for the long range cor- 
relations and for some intrinsic fluctuations in 
the hadronizing system one should rather use the 
nonextensive Tsallis statistics G|, in which one 
new parameter q describes summarily the possi- 
ble departure from the usual BG case (which is 
recovered in the q —> 1 limit) 0. Here we shall 
provide detail description of the rapidity spectra 
of particles (mostly pions) produced in pp colli- 
sions using the anomalous diffusion approach to 
this problem and in this way accounting for their 
non-equilibrium character. 

This method originates from diffusion model 
approach to nuclear multiparticle production col- 
lisions developed in jjj, which has been further 
(successfully) applied to the recent RHIC data in 
j4|. The first attempt to extend it also to the 
case of anomalous diffusion (which corresponds 
to nonextensive q ^ 1 case) has been presented 
recently for nuclear collisions || . The non- linear 
Fokker-Planck (FP) equation used in this case has 
form (/ = f(y, t) with y being rapidity and t time 



variable) : 

T^ = Ty{ J{v)F + D Ty r 



(i) 



D and J are diffusion and drift coefficients respec- 
tively The hadronization process is visualized 
here as diffusion in the rapidity space starting 
with rapidities (in cms frame) Y^h^ x ~ ± In 
(where y/s is invariant energy of reaction and 
rriT = y m 2 + (pr) 2 mean transverse mass kept 
here as given. In the approaches presented in 
H and Q linear or constant drift coefficients has 
been used and in effect obtained double gaussian- 
like form of rapidity spectra^. As was demon- 
strated in for drift proportional to the lon- 
gitudinal momentum of the particle (i.e., J(y) ~ 
sinhy) one gets the thermal (Boltzmann) distri- 
bution (for linear FP equation, i.e., for q = v = 1) 
whereas for nonlinear-FP equation the time de- 
pendent solution has the specific power-like form 
the norm of which is conserved only for \x = 1, 
therefore v = 2 — q and our distribution is given 



1 Interesingfy enough, such spectra for were already 
postulated and used on purely phenomenological grounds 
as simple parametrizations of results of string models al- 
lowing for fast numerical calculations of cosmic ray cas- 
cades 
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by the following formula: 



tot 

(1-?)— cosh(y-?/ m ) 



(2) 



2. RESULTS 



We show now results of fitting experimental 
data on rapidity distributions of secondaries pro- 
duced in pp collisions by using formula (|^) in 
the following form: 




Figure 1. Rapidity spectra obtained by UA5 
§ and Tevatron fo| fitted by formula (§). 
The broadest spectrum (full line without points) 
shows our extrapolation to the LHC energy range 
(5 TeV). 



f q (y) = f i q +) (y) + f l 

where 



(-) 



(y) 



(3) 



J a 



(y) = 



1- (!-<7)— cosh(y±y m ) 



(4) 



with T, y m and q being the (energy dependent) 
parameters (the values of transverse mass Tot was 
in this case taken as fixed and equal to tot = 0.4 
GeV, its possible energy dependence (for exam- 
ple such as used in (Q) below) is then hidden in 
the energy dependence of the parameter T) . The 
results are shown in Fig. 1 (for UA5 Q and Teva- 
tron data) where the following energy depen- 
dences of T, y m and q were used: 

T{sfs) = ai + b\ ■ In (a/s) , 



y m (Vs) = a 2 
q(yfs) = a 3 



b 2 ■ In [y/s) , 
b 3 ■ In (y/s) . 



(5) 



The values of coefficients are: ai = —0.178 ± 
0.099, bi = 0.267 ±0.016, a 2 = 0.451 ± 0.129, 
6 2 = 0.239 ± 0.023, a 3 = 0.91 ± 0.012, b 3 = 
0.072 ±0.002. It turns out that the same parame- 
ters can also fit data from P238 and UA7 @ 
groups obtained for different rapidity regions at 
energy 630 GeV, cf. Fig. 2. It is interesting that 
although our y m is not supposed to be a priori 
connected with the values of the maximal possi- 
ble rapidities at given energy, Y max , nevertheless 
we have found that ratio y m /Y max changes only 
weakly with energy varying between 53 and 1800 
GeV from 0.265 to 0.285. Notice also that for 
small values of |j/ — y m \ Eq.(||) can be aproxi- 
matcd by double quasi-gaussian form (which be- 
comes pure double gausian function in the limit 
of q^l). 




Figure 2. Comparison of data for rapidity distri- 
butions at y/s = 630 GeV obtained by P238 Jll[] 
(open circles) and UA7 Jl^] collaborations (full 
circles) with predictions of our model. 

It is well known that the charged particle mul- 
tiplicity distributions may well be fitted by a neg- 
ative binomial distribution (NB) ||. It has two 
parameters: mean charged multiplicity h, and the 
parameter k (k > 1) affecting its shape (width) 
(for k — ■> 1 NB approaches geometrical distribu- 
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tions whereas for k^ 1 — > it approaches Poisson 
distribution) . The observed widening of the nor- 
malized multiplicity distribution with increasing 
energy implies then that parameter k decrreases 
with energy. Following ideas expressed in we 
would like to bring ones attention to the fact that 
the value of parameter fc _1 may be understood 
as the measure of fluctuations of mean multiplic- 
ity. When one starts with Poisson multiplicity 
distribution and then allows n to fluctuate ac- 
cording to gamma distribution with normalized 
variance given by D(n) then it is easy to show 
|2| that as result one gets NB distribution with 
k- 1 = D{n) = a 2 (n)/(n) 2 =q-l. Namely: 
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Figure 3. The full line shows the energy depen- 
dence of obtained from NB fits to multiplicity 
distributions at collider energies ||]. The dashed 
line shows the energy variation of our parameter 
q obtained from our fits in Fig. 1. 



P(n) 



771 



r (k + n) 



T(k) 



-dn - 



r(i + n)r(fc) ( 7 + i) 



k-\-n ' 



(6) 



where 7 = k/(n). It is worth to notice that, in- 
deed, the energy variation of parameter q is al- 
most the same as fc _1 (see Fig. 3). The small 
discrepancies observed there, namely the fact that 
q — 1 is consistently larger than k~ l , can be ex- 
plained by realizing that our q has been obtained 



from fitting the experimental distributions dN / dy 
and therefore it contains in addition to the above 
mentioned fluctuations in n also fluctuations in 
inelasticity K , not discussed here. In any case it 
is at this moment tempting to assume that the 
nonextensive parameter q (parametrizing already 
nonlinearity of our FP equation used to describe 
data) also describes those fluctuations. 

In Fig. 4 we show the total inelasticity K = 
K(s) obtained by integrating spectra given by eq. 
(H) with parameters obtained from the fit to the 
corresponding rapidity distributions 



K = 



fq(y) ■ m T coshydy 



(7) 



Jo 2 

as function of the invariant energy of the reaction 
y/s both for the energies considered here and ex- 
trapolated to the LHC energy. In this case the 
values of transverse mass mx for a given energy 
was obtained by using simple interpolating for- 
mula: m T = 0.3 + 0.044 ln(Vs/20). As one can 
see in Fig. 4 such inelasticity K(s) is decreasing 
with the energy. 
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Figure 4. The energy dependence of total inelas- 
ticity obtained from our model. 

3. SUMMARY AND CONCLUSIONS 

We have provided here description of rapidity 
spectra of particles produced at CERN and Fer- 
milab energies treating their formation as diffu- 
sion process in rapidity space. To account for the 
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anomalous character of such diffusion, the non- 
linear form of Focker-Planck equation has been 
used with nonlinearity described by parameter 
q, the same as the nonextensivity parameter de- 
scribing the underlying Tsallis statistics. As seen 
in Figs. 1 and 2, very good agreement with data 
has been obtained with apparently three param- 
eters: the "temperature" T, position of the peak 
at rapidity y m and parameter q - all logaritmi- 
cally dependent on the energy of reaction ^/s [| 
However, after closer inspection it turns out that 
parameter q, which according to ||] can be re- 
garded as a measure of fluctuations existing in 
the physical system under consideration, follows 
essentiall the fluctuations of multiplicity of parti- 
cles produced at given energy. The small differ- 
ences noticed in Fig. 3 are, in our opinion, caused 
entirely by the fluctuations in the inelasticity of 
the reaction, which makes the initial energy avail- 
able for the production of secondaries a fluctuat- 
ing quantity^. Similarly, parameter y m seems to 
be closely connected with the maximal available 
rapidity Y rnax (^fs) This shows that the only 
parameter which is entirely "free" is the limiting 
temperature T 0. Therefore we can say that what 
we are proposing here is essentially one-parameter 
fit successfully describing data on rapidity distri- 
butions. 

The resultant inelasticity K (defined by the for- 
mula (Q)) turns out to be decreasing function of 
energy, as seen in Fig. 4, confirming our previ- 
ous findings in that matter [pllfl . We would like 
to bring ones attention to the fact that in our ap- 
proach we can fit equally well (using the same set 
of parameters) data obtained by P238 [jll] and 
UA7 fl2| collaborations which cover different re- 



2 We have checked that we could not get such good fits 
within formalism with q = 1, what means that nonlin- 
earity expressed by q > 1 here is essential feature of the 
data. 

3 However, we did not attempt here to estimate such fluc- 
tuations. We are planning to do it elesewhere. 
4 There is therefore possibility that it is, in fact, connected 
with Y max (K ■ \fs), i.e., it depends in the indirect way on 
the inelasticty K of the reaction. We have not pursued 
this problem here. 

5 One should notice at this point that its energy depen- 
dence as given in eq. (El) contains in itself also the possible 
energy dependence of the mean transverse mass tut, not 
accounted for in the formula (p|). 



gions in rapidity. This makes extrapolation of our 
formula to higher rapidity region more credible 
and therefore puts more weight on the obtained 
energy behaviour of inelasticity. 
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